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Abstract 



This chapter discusses an apphcation of the random matrix theory in the con- 
text of estimating the bipartite entanglement of a quantum system. We dis- 
til ' cuss how the Wishart ensemble (the earliest studied random matrix ensemble) 
appears in this quantum problem. The eigenvalues of the reduced density ma- 
trix of one of the subsystems have similar statistical properties as those of the 
Wishart matrices, except that their trace is constrained to be unity. We focus 
here on the smallest eigenvalue which serves as an important measure of en- 
tanglement between the two subsystems. In the hard edge case (when the two 
subsystems have equal sizes) one can fully characterize the probability distri- 
bution of the minimum eigenvalue for real, complex and quaternion matrices 
of all sizes. In particular, we discuss the important finite size effect due to the 
fixed trace constraint. 
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1.1 Introduction 

The different chapters of this book have already illustrated numerous appli- 
cations of random matrices in a variety of problems ranging from physics to 
finance. In this chapter, I will demonstrate yet another beautiful application 
of random matrix theory in a bipartite quantum system that is entangled. En- 
tanglement has off late become a rather fashionable subject due to its applica- 
tions in quantum information theory and quantum computation. Entanglement 
serves as a simple measure of nonclassical correlations between different parts of 
a quantum system. The more the entanglement between two parts of a system, 
better it is for the functioning of algorithms of quantum computation. This 
is so because, intuitively speaking, quantum states that are highly entangled 
contain more informations about different subparts of the composite system. 
In this chapter I will discuss the statistical properties of the entanglement in a 
particularly simple model of the 'random pure state' of a bipartite system. We 
will see how random matrices come into play in such a system. 

Indeed, historically the earliest studied ensemble of random matrices is the 
Wishart ensemble (introduced by Wishart |Wis28j in 1928 in the context of 
multivariate data analysis, much before Wigner introduced the standard Gaus- 
sian ensembles of random matrices in the physics literature) . Wishart matrices 
have found wide applications in a variety of systems (see the discussion later 
and also chapter 28 and chapter 40 of this book). In this chapter, we will see 
that the Wishart ensemble ( with a fixed trace constraint) also appears quite 
naturally as the reduced density matrix of a coupled entangled bipartite quan- 
tum system. The plan of this chapter, after a brief introduction to Wishart 
matrices, is to explore this connection more deeply with a particular focus on 
the statistics of the minimum eigenvalue which serves as a useful measure of 
entanglement. 

Let us start with a brief recollection of the Wishart matrices. Consider a 
square ( A^ x A^) matrix W of the product form W = XX^ where X is a (A^ x M) 
rectangular matrix with real, complex or quaternion entries and X' its conju- 
gate. The matrix W has a simple and natural interpretation. For example, let 
the entries Xij of the X matrix represent some data, e.g., the price of the i-th 
commodity on, say, the j-th day of observation. So, there are A^ commodities 
and for each of them we have the prices for M consecutive days, represented by 
the (A^ X M) array X. Thus for each commodity, we have M different samples. 
The product (A^ x A^) matrix W = XX' then represents the (unnormalized) 
covariance matrix, i.e., the correlation matrix between the prices of A^ commodi- 
ties. If the entries of X are independent Gassian random variables chosen from 



the joint distribution P[{Arjj}] oc exp —^Tt{XX^) 
/? = 1, 2, or 4 corresponds respectively to real, comp 



(where the Dyson index 
ex or quaternion entries), 



then the random covariance matrix W is called the Wishart matrix |Wis28| . 
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This ensemble is also referred to as the Laguerre ensemble since its spectral 
properties involve Laguerre polynomials |Bro651 IFor93| . 

As mentioned earlier, since its introduction Wishart matrices have found 
an impressive list of applications. They play an important role in statisti- 
cal data analysis |Wil621 IJohOlj , in particular in data compression techniques 
known as Prinicipal Component Analysis (PCA) (see chapter 28 and chapter 40 
of this book) with applications in image processing |Fuk90] , biological microar- 
rays [HolOOiEItOO] . population genetics |Pat06llNov08] . finance [BouOlilBurM] . 
meteorology and oceanography [Pre88] amongst others. In physics, Wishart 
matrices have appeared in multiple areas: in nuclear physics |Fyo97| , in low en- 
ergy QCD and gauge theories |Ver94a) (see also chapter 32 of this book), quan- 
tum gravity |Amb941 IAke97] and also in several problems in statistical physics. 
These include directed polymers in a disordered medium [JohOO] , nonintersect- 
ing Brownian excursions |Kat03t ISch08| and fluctuating nonintersecting inter- 
faces over a solid substrate |Nad09] . Several deformations of Wishart ensemble, 
with multiple applications, have also been studied in the literature |Ake08| . 

The Wishart matrix W has N non-negative random eigenvalues denoted 
by {wi,W2, ■ ■ ■ jWjy} (wi > for each i) whose spectral properties are well 
understood and some of them will be briefly reviewed in section 1.2. These 
include the joint distribution of N eigenvalues, the average density of eigen- 
values and also the distribution of extreme eigenvalues (the largest and the 
smallest). In this chapter we will be mostly concerned with the distribution 
of the smallest eigenvalue Wmm = min(wi, W2, • • • ,W]\f) in the particular case 
M = N corresponding to the so called hard edge (at the origin) case where the 
average (wjam) — )■ as A^ — )• oo. In this case, the properties of the small 
eigenvalues (near w = 0) are governed by Bessel functions in the large N 
limit |Ede88l IFor93t |Nag93[ IVer94b[ |Nag95| . Such hard edge properties are ab- 



sent in the traditional Wigner-Dyson Gaussian random matrix [Meh04) whose 
eigenvalues can be both positive and negative. 

The reason we are interested in the smallest eigenvalue distribution of the 
Wishart matrix is because of its application in the seemingly unrelated quantum 
entanglement problem which is the main objective of this chapter. As we will see 
later, Wishart matrices will appear naturally as the reduced density matrix in 
a coupled bipartite quantum system that is in an entangled random pure state. 
There is a slight twist though: the Wishart matrix in this system satisfies 
a constraint, namely its trace is fixed to unity. This fixed trace ensemble is 
thus analogus to the microcanonical ensemble in statistical mechanics while the 
standard (unconstrained) Wishart ensemble being the analogue of the canonical 
ensemble in statistical mechanics (for other discussions on fixed trace ensembles 
see chapter 14, section 14.3.2 of this book). In particular, our emphasis will be 
on the distribution of the smallest eigenvalue Amin in this fixed trace Wishart 
ensemble. This is because the smallest eigenvalue turns to be a very useful 
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observable in this system which contains informations about entanglement. For 
the special case M = N (hard edge), we will see that the distribution of Amin 
can be exactly computed for all N in this fixed trace ensemble in all three 
cases /3 = 1, /3 = 2 and /3 = 4. In particular, we will discuss how the fixed trace 
constraint modifies the distribution of Amin from its counterpart in the canonical 
Wishart ensemble. We will see that the global fixed trace constraint gives rise to 
rather strong finite size effects. This is relevant in the quantum context where 
the subsystems can be just a few qubits. So, it is actually important to know the 
distribution of entanglement for finite size systems (the thermodynamic limit is 
not always relevant in this context). Hence, the fact that one can compute the 
distribution of the minimum eigenvalue exactly for all N (not necessarily large) 
in presence of the fixed trace constraint becomes important and relevant. 

This rest of the chapter is organized as follows. In section 11.21 we briefly 
review some spectral properties of unconstrained Wishart matrices. In section 
11.31 we introduce the problem of the random pure state of an entangled quan- 
tum bipartite system. Its connection to Wishart matrices with a fixed trace 
constraint is established. Next we focus on the smallest eigenvalue and derive 
its probability distribution for the bipartite problem in section 11.41 Finally we 
conclude in section [LSl with a summary and open problems. 

1.2 Spectral Properties of Wishart Matrices: A brief 
summary 

Let us first briefiy recall some spectral properties of the {N x A^) Wishart ma- 
trix W = XX' with X being a rectangular (A^ x M) matrix with real (/3 = 1), 
complex (/3 = 2) or quaternion (/? = 4) Gaussian entries drawn from the joint 
distribution P[{Xij}] oc exp —^Tt{XX^) . These results will be useful for the 
problem of the random pure state of the bipartite system to be discussed in the 
next section. 

Joint distribution of eigenvalues: The N eigenvalues of W, denoted by 
{'Wi,W2, ■ ■ ■ , wn}, are non-negative and have the joint probability density func- 
tion (pdf) |Jam64] 

N 
4=1 j<k 

where a = {\+M — N) — 2/ 13 and the normalization constant K^^m can be com- 
puted exactly | Jam64j . Without any loss of generality, we will assume A^ < M. 
This is because if A'^ > M, one can show that N — M eigenvalues are exactly 
and the rest of the M eigenvalues are distributed exactly as in Eq. ()1.2.ip 
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with A^ and M exchanged. Note that while for Wishart matrices M — A^ is a 
non-negative integer and /3 = 1, 2 or 4, the joint density in Eq. (jl.2.ip is well 
defined for any (3 > and a > — 2//3 (this last condition is necessary so that 
the joint pdf is normalizable) . When these parameters take continuous values 
the joint pdf is called the Laguerre ensemble. 

Coulomb gas interpretation, typical scaling and average density of 
states: The joint pdf (|1.2.ip can be written in the standard Boltzmann form, 
P[{u;i}] oc exp [—(3E{{wi})] where 

1 ^ 1 

E[{wi} = 2 X]("^* -alogWi) - -^ln\wj - Wk\ (1.2.2) 

can be identified as the energy of a Coulomb gas of charges with positions 
{wi} > 0. These charges repel each other via the 2-d Coulomb (logarithmic) 
interaction (the second term in the energy), though they are restricted to live 
on the positive real line. In addition, these charges are subjected to an external 
potential which is linear+logarithmic (the first term in the energy ). The exter- 
nal potential tends to push the charges towards the origin while the Coulomb 
repulsion tends to spread them apart. The first term typically scales as WtypN 
where Wtyp is the typical value of an eigenvalue, while the second term scales 
as N'^ for large N. Balancing these two terms one gets u^typ ~ ^ for large A^. 
Indeed, this scaling shows up in the average density of states (average charge 
density) which can be computed from the joint pdf and has the following scaling 
for large N 

1 ^ 1 

Pn{w) = - ^{6iw - Wi)) ^ - /MP (^) (1.2.3) 

where the Marcenko-Pastur(MP) scaling function is given by |Mar67| (see also 
chapter 28 section 28.4.1 of this book) 

1 



/mp(x) = —- ^/{h-x){x-a). (1.2.4) 

zvrx 

Thus the charge density is confined over a finite support [a, b] with the lower 
edge a = (1 — c^^/^)^ and the upper edge b = {1 + c~^l'^)'^ with < c = 
N/M < 1. For all c < 1, the average density vanishes at both edges of the 
MP sea. For the special case c = 1 (this happens in the large A^ limit when 
M — N « 0{N)), the lower edge a gets pushed towards the hard wall at 
(this is the so called hard edge limit) and the upper edge 6—7-4 and the average 
density simply becomes, fupix) = ^^/i^ — x)/x for < x < 4. It diverges as 
x~^''^ at the hard lower edge x = 0. 

For later purposes, it is also useful to calculate the average value of the trace 
Tr = Ylii=i "^i- Using the expression for the average density of states, it follows 
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that for large N 

lTr)=N wpN(w)dw-^ — . (1.2.5) 

Jo c 

In particular, for c = 1 (i.e., M — N « 0{N)), we have (Tr) = N"^ in the 
large N limit. 

Maximum eigenvalue: Let t^max = max{wi,W2, ■ ■ ■ ,wn) denote the maxi- 
mum eigenvalue. On an average, it is located at the upper edge of the MP den- 
sity of states. It then follows from Eq. (|1.2.3p that (tfrnax) = bN for large A^. 
However, for large but finite A^, the random variable Wmax fluctuates, from one 
sample to another, around its average value bN. The typical fluctuations around 
its mean were shown to be ~ 0{N^'^) for large N |Joh001 IJoh01| and the lim- 
iting distribution of these typical fluctuations are described by the well known 
Tracy- Widom density |Tra94) . In other words, Wmax = bN + c^l^b'^/^N^/'^Xi 
where the random variable x has an A^- independent limiting pdf, 5/3 (x) de- 
scribed by the Tracy- Widom function |Tra94| . In contrast, atypically large, 
e.g., ~ 0{N) fluctuations of t^max from its mean are not described by the TW 
density. Such large fluctuations play an important role in many practical ap- 
plications such as in PCA ^Viv07t Maj09a| . Far away from the mean 6A^, these 



atypically large fluctuations of Wmax are instead described by large deviation 
functions associated with the pdf of P{wyaa.x-, A^) and are of the form 



P{wma.x = t,N) ~ exp 
~ exp 



t-bN 



-/3Af$_ 



A^ 



for t « bN (1.2.6) 
for t » bN. (1.2.7) 



The left rate function ^^(x) was computed explicitly for all c in |Viv07| ex- 
tending a Coulomb gas approach developed originally in |Dea06| to compute 
the corresponding left rate functions for Wigner-Dyson Gaussian matrices. On 
the other hand, the computation of the right rate function $+(x) required a 
different approach and was recently obtained explicitly for all c |Maj09a] . The 
right rate function in the Wigner-Dyson Gaussian case was also obtained by 
a different, albeit rigorous, method in |Ben01| . One interesting point is that 
while the limiting TW density g/six) depends on /3, the rate functions ^zp{x) 
are independent of /3. 

Minimum eigenvalue: Since in this chapter our main interest in the problem 
of bipartitite entanglement concerns the lowest eigenvalue of the reduced density 
matrix, we need to discuss, in some detail, the statistical properties of the 
minimum eigenvalue of the unconstrained Wishart ensemble. For the minimum 
eigenvalue, u'min = rain{'Wi,W2, . . . ,wn), the situation is rather different for 
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c < 1 and c = 1 cases. For c < 1, the lower edge of the MP sea is at a = 
(1 — c^"*^")^ > 0, indicating that (wmm) = o.N in the large A^ limit and thus 
the typical value of Wram ~ 0{N). The typical fluctuations of Wmin around 
this mean value aN are again described by the TW density (appropriately 
rescaled). The large deviation functions describing atypical fluctuations, to our 
knowledge, have not been systematically studied as in the maximum eigenvalue 
case (though see |Che96| and references therein). 

The situation, however, is quite different in the c = 1 case (when M — 
N « 0{N) where the lower edge of the MP sea a — >• 0. This is the so called 
hard edge case. We will see shortly that in this case the typical value of the 
minimum eigenvalue scales as u^min ~ 1/-^ for large N, to be contrasted with 
the behavior Wuiin ~ aN for c < 1. There have been a lot of studies on the 
distribution of Wmin in this hard edge c = 1 (M — N « 0{N)) case, notably 
by Edelman [Ede88j and Forrester [For 931 IFor94] . It has also found very nice 
applications in QCD (see e.g. chapter 32, section 32.2.6 of this book). Here, for 
simplicity, we will focus on the special case M = N (such that c = 1 strictly for 
all N, and not just for large A''). For other cases when M — N ^ 0(1) (so that 
c = 1 only in the large N limit), a summary can be found in the table 32.2 of 
chapter 32 of this book (see also section 4.2 of |Ake08] and references therein). 

In this special case M = N, the cumulative distribution of the minimum, 
Qn{z) = Prob[^i;min > z,N], is known |Ede88] exactly for all N in all the three 
cases /? = 1, /3 = 2 and /3 = 4. Note that, Qn{z) = J^ ■ ■ ■ J^ P[{wi}] U dwi 
where ^[{TUi}] is the joint pdf given in Eq. (11.2.ip . For M = N, this multiple 
integral Qn{z) can be easily performed for /3 = 2 by making a trivial shift 
Wi ^ Wi + z and one gets for all N 

QAr(z) = exp[-iVz]; /3 = 2 (1.2.8) 

For /3 = 1 and 4, the simple shift does not work. However, the integral Qn{z) 
can be calculated explicitly [Ede88] . For /3 = 1, one obtains for all A'' 

Q,(,) . /(^ + i> f r^n ,-»./. u(^i^.-'-A,y (1.2.9, 

^^' 2^-y^T{N/2) J^ V 2 ' 2'2y ^ ^ ^ 

where U{p,q,z) is the confluent (Kummer) hypergeometric function |Abr72] . 
For /3 = 4, while Edelman does not provide an explicit expression for Qm{z), it 
is not difficult to obtain Qn{z) by using his Lemma 9.2 [Ede88] and one gets 
(see also |For94] ) 

Qiv(^)=e-2^^Fi('-Ar;i;-z'); /3 = 4 (1.2.10) 

where 

,F,(p;,;.) = l + ?£ + ?fc±i| + ... (1.2.11) 

q 1! q{q+ 1) 2! 
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is the degenerate hypergeometric function |Abr72] . 

The large A^ limit is interesting where in all three cases /3 = 1, 2 and 4, 
the cumulative distribution of the minimum Qn{z) approaches a scaling form: 
Qn{z) — )• q/3{zN) where the scaling function ^/^(y) can be computed explicitly 



<ii{y) 
52 (y) 

54 (y) 



exp[-^-y/2] 
exp[-y] 

1 ^ ^-2y+2^ _^ ^-2y-2^' 



Note in particular that the average value {w^ 
large A^, where the prefactor 



!^ QN{z)dz 



C/3 



qp{y)dy 



(1.2.12) 
(1.2.13) 

(1.2.14) 
cp/N for 

(1.2.15) 

(1.2.16) 
(1.2.17) 
(1.2.18) 



where erf (z) = —i= Jq e " dn is the standard error function and erfc(2;) 
1 — erf (2:). These results will be used in Section [T7S1 



can be computed explicitly in all three cases and one gets 



Cl = 


= 2 


1 - y^erfc(l/^/2) 


= 0.68864.. 


C2 = 


= 1 




C4 = 


1 

2 


1 + y^erf(l/^/2) 


= 1.20534.. 



1.3 Entangled Random Pure State of a Bipartite 
System 

We now turn to the main problem of interest in this chapter, namely the prop- 
erties of an entangled random state of a quantum bipartite system. We will see 
that Wishart matrices, albeit with a fixed trace constraint, play a central role 
in this problem. 

As mentioned in the introduction, entanglement has been studied exten- 
sively in the recent past due to its central role in quantum information and 
possible involvement in quantum computation. In the context of quantum al- 
gorithms, it is often desirable to create states of large entanglement. A potential 
candidate for such a state with 'large entanglement' that is relatively simple to 
analyse turns out to be the 'random pure state' in a bipartite system |Hay06| 
which we will describe in detail shortly. Such a random pure state can also 
be used as a null model or reference point to which the entanglement of an 
arbitrary time-evolving state may be compared. Apart from the issue of bi- 
partite entanglement, statistical properties of such random states are relevant 
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for quantum chaotic or non-integrable systems. The apphcabihty of random 
matrix theory and hence of random states to systems with well-defined chaotic 
classical limits was pointed out long back |Boh84| . 

We start with a general discussion of entanglement in a bipartite setting 
without any reference to any specific statistical measure. The statistical prop- 
erties will be discussed later when we introduce the 'random' state. For now, 
the discussion below holds for any quantum pure state. Let us consider a com- 
posite bipartite system A0 B composed of two smaller subsystems A and B, 
whose respective Hilbert spaces "H^ and T-L^ have dimensions N and M. 
The Hilbert space of the composite system 7^(^*^) = Ji^^ T-L^ is thus 
ATM-dimensional. Without loss of generality we will assume that N < M. 
Let {|^ )} and {|a^)} represent two complete basis states for A and B respec- 
tively. Then, any arbitrary pure state IV') of the composite system can be most 
generally written as a linear combination 

N M 



5]5]x,,„|i^)(»|a^) (1.3.1) 



j=l a=l 

where the coefficients Xj^q's form the entries of a rectangular {N x M) matrix 
X = [xi^a]. As an example of such a bipartite system, A may be considered 
a given subsystem (say a set of spins) and B may represent the environment 
(e.g., a heat bath). 

Next we discuss the density matrix and the concept of entanglement. For a 
pure state, the density matrix of the composite system is simply defined as 

P=|V')(^I (1.3.2) 

with the constraint Tr[/3] = 1, or equivalently (V'l^) = 1- Note that had the 
composite system been in a statistically mixed state, its density matrix would 
have been of the form 

P = X]^'^I^^^^^'=I' (1.3.3) 

k 

where |V'fc)'s are the pure states of the composite system and < pfc < 1 
denotes the probability that the composite system is in the fc-th pure state, 
with '^j.Pk = 1. A classical example of such a mixed state is when the system 
is in the canonical ensemble at given temperature T: in this case the density 
matrix is given by 

P = Y.^e-^"''^^\E){E\ (1.3.4) 

E 

where Z = ^^ ^-E/kBT jg ^j^g canonical partition function (/c^ is the Boltz- 
mann constant) and the pure state E) denotes the energy eigenstate (with 
eigenvalue E) of the full system. In this chapter, we will not discuss the mixed 
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state and will restrict ourselves only to the case of a pure state whose density 
matrix is given in Eq. (jl.3.2p . 

The concept of entanglement is simple. A pure state l'^) is called entangled 
if it is not expressible as a direct product of two states belonging to the two sub- 
systems A and B. Only in the special case when the coefficients have the prod- 
uct form, Xi^a = o.iba for all i and a, the state \ijj) = \(l) ) ® I'P^) can be written 
as a direct product of two states \4> ) = X^j=i ai\i ) and \(j)^) = ^^=1 ba\a^) 
belonging respectively to the two subsystems A and B. In this case, the com- 
posite state \ijj) is fully separable or unentangled. But otherwise, it is generically 
entangled. 

Upon using the decomposition in Eq. (|1.3.ip . the density matrix of the pure 
state can be expressed as 

'° = EE^*."4/3l^^)0'^l»l«'')(/3''l (1.3.5) 

where the Roman indices i and j run from 1 to A^ and the Greek indices a and (3 
run from 1 to M. We also assume that the pure state \tp) is normalized to unity 
so that Tv[p\ = 1. Hence the coefficients Xj_q,'s must satisfy Ylia kj,«P = 1- 

Given the density matrix of the pure composite state in Eq. (J1.3.5P , one can 
then compute the reduced density matrix of, say, the subsystem A by tracing 
over the states of the subsystem B 

M 

PA = TTB[p\ = Y,{a''\p\a''). (1.3.6) 

a=l 

The reduced density matrix is important because if we measure any observable 
O of the subsystem A, its expected value is given by Tr[Opyi]. Thus, pA is the 
basic physical object whose properties are directly related to measurements. 
Using (J1.3.5P one gets 

N M N 

P^ = E E^^."4«N^)0-^l = E W^,\^''){A (1.3.7) 

i,j=la=l i,j=^ 

where Wij 's are the entries of the N x N square matrix W = XX' . In a similar 
way, one can express the reduced density matrix ps = Tia [p] of the subsystem 
B in terms of the square M x M dimensional matrix W = X^X. Hence we 
see how the Wishart covariance matrix W = XX' appears in this quantum 
problem. 

Let Ai, A2, . . . , Atv denote the N eigenvalues oi W = XX^ . Note that these 
eigenvalues are non-negative, Aj > for all i = 1,2,... , A. Now the matrix 
W = X^X has M > N eigenvalues. It is easy to prove that M — A of them 
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are identically and N nonzero eigenvalues of W are the same as those of W. 
Thus, in this diagonal representation, one can express pA_ as 

N 

PA = Y.X^\Xt){Xt\ (1.3.8) 

where |A^)'s are the eigenvectors of T^ = XX' . A similar representation holds 
for pb- It then follows that one can represent the original composite state {ip) 
in this diagonal representation as 

N 



where |A^) and \Xf) represent the normalized eigenvectors (corresponding to 
the same nonzero eigenvalue Aj) ofW = XX' and W = X'X respectively. This 
spectral decomposition in Eq. (jl.3.9p is known as the Schimdt decomposition. 
The normalization condition (tplip) = 1, or equivalently Tr[p] = 1, thus imposes 
a constraint on the eigenvalues, X^j^^ Ai = 1. 

Note that while each individual state |A|^) ® |Af ) in the Schimdt decom- 
position in Eq. ()1.3.9p is separable, their linear combination {tp), in general, is 
entangled. This simply means that the composite state \ip) can not, in general, 
be written as a direct product \ip) = |(/)^) \(l)^) of two states of the respective 
subsystems. The spectral properties of the matrix W, i.e., the knowledge of 
the eigenvalues Ai, A2, . . . , Xn, in association with the Schimdt decomposition 
in Eq. (|1.3.9p . provide useful information about how entangled a pure state is. 

Measures of Entanglement: It is useful to construct a measure of entan- 
glement, i.e., a function of the eigenvalues Aj's whose value will tell us how 
entangled a pure state is. There are many ways of constructing such a mea- 
sure. Its value should monotonically increase from the configuration of Aj's 
where the state is fully separable to the configuration where the state is max- 
imally entangled. These two configurations, recalling that Yli Xi = 1, are the 
following: 

(i) separable: When one of the eigenvalues, say Ai is 1 and the rest are all 
identically zero. Then the state completely decouples as only one term, say the 
first term, is present in Eq. ()1.3.9p . 

(ii) maximally entangled: When all eigenvalues are equal, i.e., Xi = !/N. In 
this case all A^ terms in Eq. (J1.3.9P are present. 

In Fig. (jl.ip . we present a cartoon for A^ = 3 for the purpose of illustra- 
tion. In the three dimensional space (Ai,A2,A3), any point on the triangular 
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plane Ai + A2 + A3 = 1 with Aj > represents an allowed configuration. The 
three vertices, where the system gets completely factorised, represent the fully 
separable configurations (situation (i) above). On the other hand, the centroid 
(1/3, 1/3, 1/3) represents the fully (maximally) entangled configuration (situa- 
tion (ii) above). 



sepaS^ 



maximally entangled 




A. 1/(1, 0,0) 
separable 



Figure 1.1: A cartoon for A^ = 3. The system lives on the triangular plane Ai + 
A2 + A3 = 1. The vertices of the triangle represent the separable configurations 
and the centroid (Ai = A2 = A3 = 1/3) represents the maximally entangled 
configuration. 



^Ar 



If the system is in a given configuration {Aj} on the allowed plane ^i=i Aj = 
1 (and Aj > 0), how much entangled it is? In other words, how do we mea- 
sure the entanglement content of a given configuration of {Aj}? This is usually 
done by defining the so called entanglement entropy, a single scalar quantity 
associated with each configuration, i.e., each point on the plane X^j=i Aj = 1 
(and Aj > 0). A standard and perhaps most studied measure of entangle- 
ment is the so called von Neumann entropy |Ben96] . Si = — X]j=i Ajln(Aj) = 
— Tr[/5^1n(/9A)], which has its smallest value 5i = in configuration (i) and 
its maximum possible value Si = InA^ in configuration (ii). Renyi entropy de- 
fined as Sq = In (^^ A^) /(I — q) |Ren70] with the parameter q' > is a natural 
generalization that reduces to the von Neumann entropy when g — )■ 1. Again, 
for any q > 0, Sq = at the 'separable' vertices (situation (i) above) and 
Sq = In(A^) at the 'maximally entangled' centroid (situation (ii) above). For 
q = 2, X]j=i '^j = exp[— ^2] is called the purity that has been widely studied 
(see [FacOG] and references therein). For other measures we refer the reader 
to the introduction in j GirOTj . Essentially, one can define any scalar quantity 
whose value increases monotonically as one moves from fully 'separable' to max- 
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imally 'entangled' configurations (e.g., as one moves from tlie vertices towards 
tlie centroid in Fig. (jl.ip ) (for more detailed prescriptions and requirements on 
the measure, see e.g. ||Ved98j ). 

Important informations regarding the degree of entanglement can also be ob- 
tained from the two extreme eigenvalues, the largest Amax = max(Ai, A2, . . . , A^v) 
and the smallest Amin = min(Ai, A2, • • • , Xn)- Due to the constraint J2i=i ^j = 1 
and the fact that eigenvalues are all non-negative, it follows that 1/A^ < Amax ^ 
1 and < Amin < 1/A^. Consider, for instance, the following limiting situations. 
Suppose that the largest eigenvalue Amax = max(Ai, A2, . . . , Xn) takes its max- 
imum allowed value 1. Then due to the constraint Yli=i Xi = I and the fact 
that Aj > for all i, it follows that all the rest {N — 1) eigenvalues must be 
identically 0. Thus it corresponds to the configuration (i) above of fully sepa- 
rable state. On the other hand, if Amax = 1/-^ (i-e., it takes its lowest allowed 
value), it follows that all the eigenvalues must have the same value, Aj = 1/A^ 
for all i, again due to the constraint Yli=i^i — 1- This then corresponds to 
situation (ii) of maximally entangled state. Thus, for instance, one can consider 
— In(Amax) as a measure of entanglement as it increases from its value in the 
separable state to its maximal value In(A^) in the maximally entangled case. In 
fact, — In(Amax) is precisely the 5 — )• 00 limit of the Renyi entropy Sq. 

In this chapter our particular interest is on the smallest eigenvalue < 
Amin < 1/A^. When Amin takes its maximal allowed value 1/A^, it follows again, 
from the constraint X]i=i ^i = 1 and Aj > for all i, that all the eigenval- 
ues must have the same value Aj = 1/A^. This will thus make the state \ip) 
maximally entangled, i.e., situation (ii). In the opposite case, when Amin takes 
its smallest allowed value 0, while it does not provide any information on the 
entanglement of the state \ip), one sees from the Schmidt decomposition (I1.3.9P 
that the dimension of the effective Hilbert space of the subsystem A reduces 
from A^ to A^— 1 (assuming that Amin is non-degenerate). Indeed, if Amin is very 
close to zero, one can effectively ignore the term containing Amin in Eq. (I1.3.9P 
and achieve a reduced Hilbert space, a process called 'dimensional reduction' 
that is often used in the compression of large data structures in computer vi- 
sion [Wil621 IFukDOl IViv07| . Thus the knowledge of Amin and in particular its 
proximity to its upper and lower limits provide informations on both the degree 
of entanglement as well as on the efficiency of the dimensional reduction process. 

Random Pure State: So far, our discussion was valid for an arbitrary pure 
state in Eq. ()1.3.ip with any fixed coefficient matrix X = [xi^a]- One can now 
introduce a statistical measure or distribution for the entries of the matrix X 
which, in turn, will induce a probability distribution for the eigenvalues Aj's 
of W = XX'^ that appear in the Schimdt representation in Eq. (jl.S.Op . As 
a result, any measure of entanglement (e.g., the von Neumann entropy or the 
minimum eigenvalue Amin) will also have a statistical distribution associated 
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with it. The main challenge then is to compute this probability distribution of 
the entanglement, given the measure on the entries of X. 

So, what is an appropriate measure on X? Evidently, we can not choose any 
arbitrary measure on X. Indeed, in Eq. (I1.3.ip we can just rotate (unitarily) the 
basis of the Hilbert space. Clearly physical properties of the system in a pure 
state should not depend on which basis we choose. Thus the joint probability 
distribution of the entries Xi^a in Eq. (11.3. ip should be invariant under a unitary 
(or orthogonal if we restrict X to be real) transformation \ip) — )• Ulip), where U 
represents a unitary operator. The only measure that remains invariant under 
a unitary rotation is the uniform measure over all pure states. This is the 
so called Haar measure where the coefficients Xi^^s are uniformly distributed 
over all possible values satisfying the constraint ^^ ^ Ixj^qP = 1 or equivalently, 
P[{Xjj}] oc 6 (Tr[XX^ — l). The physical meaning of this Haar measure is 
clear: under unitary time evolution, and in absence of any other conservation 
law (such as fixed energy as in the case of standard microcanonical ensemble in 
statistical physics), the system visits all allowed normalized pure states equally 
likely, i.e., Haar measure is a natural stationary measure under unitary evolution 
when ergodicity holds over all allowed pure states of the composite system. This 
is, in fact, the case in many physical situations when the system is described 
by a sufficiently complex 'time-dependent' Hamiltonian as in quantum chaotic 
systems |Ban02| . 

Given that the entries of X are distributed via the Haar measure P[{Xjj}] oc 
5 (Tr(XX^^ — l) , the next question is how are the eigenvalues oi W = XX^ 
distributed? Noticing that the eigenvalues Aj's of W = XX^ are the same as 
Wishart eigenvalues, except with the additional constraint Tr(Ty) = J2i=i -^« — 
1, it follows immediately from Eq. ()1.2.ip that the joint pdf of Aj's is given 
by [LI088I IZycOT] 



/ N \ ^ B 

P[{AJ] = Bm,nS E ^^ - 1 n X!^''-''^'^-' n 1^.- - ^'^1' (1-3-10) 

\i=l / i=l j<k 

where the normalization constant Bm,n is known explicitly [ZycOl] . Note that 
the exponential factor e~2 5I,i=i-*^» present in Eq. (11.2. ip becomes a constant 
due to the constraint Y2i=i ^i = ^ ^-^d hence is absorbed in the normalization 
constant Bm,n- The ensemble described in Eq. (11.3. lOp can thus be seen as 
the microcanonical version of the canonical Wishart ensemble in (11.2. ip . 

When the coefficient matrix X is drawn from the Haar measure, we will refer 
to the state in Eq. ()1.3.ip as a random pure state. Given that Aj's corresponding 
to the random pure state are distributed via the joint pdf ()1.3.10p . it follows that 
the associated observables such as the von Neumann entropy Si = — Ajln(Aj), 
the maximum eigenvalue Amaxi the minimum eigenvalue Amin etc. are also 
random variables. The main technical problem then is to evaluate the statistical 
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properties (such as the mean, variance or even the full probability distribution) 
of such observables. 

There have been quite a few studies in this direction. For example, the 
average entropy (^i) (where the average is performed with the measure in Eq. 
(jl.3.10p ) was computed for /3 = 2 by Page |Pag95| and was found to be (5) ss 
In(A^) — 257 for large 1 << A^ < M. Noting that In(A^) is the maximal possible 
value of entropy of the subsystem A, the average entanglement entropy of a 
random pure state was concluded to be near maximal. Later, the same result 
was shown to hold for the /? = 1 case [Ban02] . On the other hand, there have 
been only few studies on the full probability distribution of the entanglement 
entropy. The distribution of the so called G-concurrence |Cap06| , a measure 
of entanglement, was computed exactly in the large A^ limit and was shown to 
have a point measure (delta function). For small N, the distribution of purity 
is known [Gir07] . On the other hand, for large N, the Laplace transform of the 
distribution of purity (for positive Laplace variable) was computed in |Fac081 
IPas09| which only gave partial information about the full purity distribution. 

Recently, using a Coulomb gas approach, the full probability distribution of 
the Renyi and von Neumann entropy, as well as that of purity, was computed 
exactly in the large N limit by studying the associated Coulomb gas model via 
a saddle point method [NadlO] . Interestingly, the pdf of the entropy exhibits 
two singular points which correspond to two interesting phase transitions in 
the Coulomb gas problem jNadlOl IFacOSl IPas09] . Similar phase transitions in 
the Coulomb gas picture, leading to a nonsingular pdf of a physical observable, 
have also been noted recently in several other problems where the random 
matrix theory is applicable: these include the pdf of the conductance and the 
shot noise power through a mesoscopic cavity |Viv08[ IVivlOl IOsi08| (see also 
the chapter 35 and 36 of this book for applications of the RMT to quantum 
transport properties), the pdf of the number of positive eigenvalues (the so 
called index) of Gaussian random matrices |Maj09b| , nonintersecting Brownian 
interfaces near a hard wall |Nad09] and in information and communication 
systems [Kaz09] to name a few. 

Here our focus is on the statistical properties of the minimum eigenvalue 
Amin and for all values of M = A^. For the special case /3 = 2 and M = N, the 
average value (Amin) was studied recently by Znidaric |Zni07j . He computed, 
by hand, (Amin) for small values of A^ and conjectured that (Amin) = l/N^ for 
all A^. Later, in |Maj08| , the full probability distribution of Amin was computed 
explicity for all M = A^ and /? = 1 and (3 = 2. Znidaric's conjecture for 
(3 = 2 then followed as a simple corollary |Maj08| . In the next section, I briefly 
outline this derivation and also provide a new result for the distribution of Amin 
for /3 = 4. 
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1.4 Minimum Eigenvalue Distribution for M = N 

In this section, we compute the distribution of Amin when the eigenvalues are 
distributed via Eq. ()1.3.10p . It is easier to compute the cumulative distribution 

i?7v(a^) = Prob [Amin > 2;] = Prob [Ai > x, A2 > x, . . . , A^r > x] . (1-4.1) 

Using (I1.3.1UD 



f(M-A.+l)-l^^^_ 



poo poo / N \ N 

Rn{x) = Bm,n •••/ -^ 5^Ai-l mi^i-^fclTI^ 

(1.4.2) 
The challenge is to evaluate this multiple integral. 
We proceed by introducing an auxiliary integral 



poo poo / N \ AT 

I{x,t)= / ••• / 6i^X,-t]l[\Xj-\kfllX, 
■^^ ■^^ \i=i / j<k i=i 



f(M-A.+l)-l 



(1.4.3) 
If we can evaluate I{x, t) for all t, then 

RNix) = BM,Nlix,l). (1.4.4) 

To evaluate I{x,t), it is natural to consider its Laplace transform 

N 



•^0 Jx Jx ■ ,,, •_i 



j<k i=l 

(1.4.5) 
Next, a change of variable Xi = -^Wi reduces it to 

/ I{x,t)e-^'dt=(l] / .../ e~^^l^-^ll\w,-w,fll 

Jo V^^y J2SX/I3 J2SX/I3 f^f. fJl 

(1.4.6) 
where a = (1 + M — N) — 2/(3. Next we recognize the multiple integral, up to 
an overall constant, as the cumulative distribution Qn{2sx/ (3) of the minimum 
eigenvalue u;min in the unconstrained Wishart ensemble discussed previously. 
Thus, up to an overall constant Ai independent of s, we have 

I" /(x, t)e~''dt = A^s-P^'^/^ Qn (^) ■ (1.4.7) 

The program then is to invert this Laplace transform, compute /(x, t) for all 
t and calculate Rn{x) using Eq. (jl.4.4p . Henceforth, we will drop the overall 
constant which can be finally fixed from the normalization that i?Ar(0) = 1. 



Wj dwi 
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Thus, if we know the cumulative distribution of the Wishart minimum eigen- 
value Qm{z), we can, at least in principle, determine the minimum eigenvalue 
distribution Rn{x) for the random pure state problem. This is hardly sur- 
prising given the microcanonical to canonical correspondence between the two 
ensembles. In practice, however, it is nontrivial to invert the Laplace transform 
in Eq. ()1.4.7p . That indeed is the real challenge. We will see below that fortu- 
nately for M = N, where QNiz) is given in Eqs. (fTXQl) . (fTXHI) and (II. 2. 101) 
for /3 = 1, 2 and 4 respectively, this Laplace inversion can be carried out in 
closed form and one can compute RNix) explicitly in all three cases ^ = 1, 2 
and 4. 

The case /3 = 2: Let us start with the simplest case /3 = 2 with M = N. Here 
Qn{z) = e~^^ from Eq. ()1.2.8p . Dropping the overall constant, Eq. ()1.4.7p 
gives 

oo —sNx 

I(x,t)e-''dt = ^^. (1.4.8) 

s'^ 

The Laplace inversion is trivial upon using the convolution theorem giving (up 

to an overall constant) I{x,t) = {t — Nx) ~^@{t — Nx) where @{z) is the step 

function. Putting t = 1 and using (|1.4.4p gives the exact distribution |Maj08| 

Rn{x) = Prob [Amin > x] = (1 - Nxf^'^ e (1 - Nx) . (1.4.9) 

Subsequently, the pdf is given by 

P^(x) = _^%M = Ar(Ar2 _ i)(i _ Nxf-^ G(l - Nx). (1.4.10) 

ax 



A plot of this pdf can be found in Fig. (jl.2p for N = 2. Thus Pn{x) in 
X G [0, 1/A^] has the limiting behavior 



Pn{x) -^ NiN"^ -I) as x^O 
-^ N{N^-1){1-Nx)^'-^ 

One can easily compute all the moments explicitly 



iV (iV^ - 1) (1 - TVx)^'-^ as x^l/N (1.4.11) 



f^kiN) = {XiJ = j^ x^Pn{x) dx = ^t^]^)^!/t)^ • (1-4-12) 

In particular, for k = 1, we obtain for all A^ 

m(iV) = (Amin) = ^, (1.4.13) 

proving the conjecture by Znidaric |Zni07| . Putting k = 2 in Eq. (|1.4.12p . we 

2 

Af4(Ar2+i)- 



get the second moment /U2 = ivttwttt- Thus the variance is given by 



1 / /v2 _ r 

a' = t.,iN) - [^,iN)f = ^ f ^^^^ ) . (1.4.14) 
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Figure 1.2: The pdf of Amin for iV = 2 for /3 = 1, /3 = 2 and /3 = 4. 



The case /? = 1: The computation in this case proceeds as in the case /3 = 2, 
though the Laplace inversion is nontrivial. Omitting details |Maj08| , we just 
quote here the main results. For M = N, the pdf of the minimum eigenvalue, 
Pn{x) = —dRN{x)/dx is nonzero in x G [0, 1/A^] and is given by |Maj08| 



Pn[x) = Am x-^/2 (i_^^)(^^+^-4)/2 ^p^ f 



N + 2 N -I N^ + N -2 I- Nx 



2 ' 2 ' 2 ' X 

(1.4.15) 

where 2-^1 (a, ^c,c, z) is the standard hypergeometric function |Abr72) and the 
constant An is given by 



A 



N 



NT{N)r{N'^/2) 



2N-1 Y(N/2) r((7V2 + AT _ 2)/2) ■ 
The limiting behavior of P/v(2;) as x — )• and x — )■ 1/A^ can be worked out 

/?r(iv)r(7vV2) 



1.4.16) 



^N[x) ^ [2^^-ir2(iv/2)r((iv-i)/2)J "^ ""' "^^^ 

« ^^iV-^/2(l-Arx)(^'+^-4)/2 ^ x^i/N (1.4.17) 



All moments can also be worked out explicitly |Maj08] . In particular, the 
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average is given by 

/ii(iV) = (Amin) = ^rr(Ar/2)r((iV + 5)/2)2^-i ^^' [^' 2' ~1~' ^ " ^ 

(1.4.18) 
Thus the expression for (Amin) for arbitrary A^ in the real (/3 = 1) case is consid- 
erably more complicated than its counterpart in Eq. (J1.4.13P for the complex 
case. One finds, from Eq. (jl.4.18p . that f^i{N) decreases with increasing A^, 
e.g.,/ii(l) = 1, /xi(2) = (4-7r)/8, /ii(3) = (2-\/3)/9 etc. One can show |Maj08| 
that asymptotically for large A^, fJ-i^N) decays as 

f^iiN) « ^3 (1.4.19) 

where ci is precisely the constant in Eq. ()1.2.16p . 

The case /3 = 4: For /3 = 4, we first substitute Qn{z) from (|1.2.10p in 
(ll.4.7p . expand the hyper geometric function in power series as in (II. 2. lip and 
then invert the Laplace transform term by term to get a series for I{x,t). To 
transform each term, we make use of the convolution theorem. We then put 
i = 1 in the expression for I{x,t) and compute i?Af(x) in ()1.4.4p . The overall 
constant is fixed by imposing -R7v(0) = 1. This gives the explicit expression, 
vahd for < X < l/N, 

The pdf PAr(x) = —dRN{x)/dx of Amin vanishes linearly at x = as Pn{x) — )■ 
Bnx where B^ = N{2N + 1)(2A^2 _ i)(2iV2 - 2)/3. This in contrast to the 
/3 = 1 case (where Pn{x) diverges as x~"^'^ as x — >• 0) and also to the (3 = 2 case 
where Pn{x) approaches a constant as x — )■ (see Fig. (|1.2p ). At the upper 
edge, when x — )■ 1/A^, the pdf vanishes as Pn{x) ~ (1 — Nx)"^^ ~N-2 ^ 

All the moments can also be calculated explicitly for /3 = 4. For example, 
the average is given by 

«(A')MA.™>^^|:(:)i;(i)'- ("-^D 

One can extract the large A^ asymptotics of this sum. We first express ( ;. ) = 
T{N+l)/V{N—k+l), then use the property of the Gamma function, lim^^oo r(z+ 
a)/T[z) — )• z", to obtain for large A^ 

1 ^ A-l 
fc=0 ^ '' 
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The sum can be exactly evaluated giving, 

/ii(iV) = (A,ni„) « ^ (1.4.23) 

where C4 is precisely the constant in Eq. (jl.2.18p . 

Let us then summarize the behavior of the pdf Pn (Amm) of Amin in the three 
cases /3 = 1, /3 = 2 and /? = 4 (see Fig. (fL2]l ). At the lower edge x -)■ 0, Pn{x) 
displays very different behavior in the three cases. As x — )• 0, PAr(x) diverges 
as x~^''^ for /? = 1, approaches a constant for /3 = 2 and vanishes linearly for 
/3 = 4. On the other hand, at the upper edge x — )• 1/A^, Pn{x) approaches as 
a power law in all three cases, albeit with different powers, Pn{x) ~ (1 — Nx^i^ 
where ui = {N"^ + N - 4)/2, U2 = N"^ - 2 and 2^4 = 2iV2 _ ^r _ 2. 

In the large A^ limit and in the range x << 1/A^ (far away from the upper 
edge) the cumulative distribution Rn{x) = jj P]\[{x')dx' approaches the 
scaling form Rn{x) — )• qp{xN'^), where the scaling functions qpiy) are exactly 
same as in the unconstrained Wishart case given respectively in Eqs. ()1.2.12p . 
(I1.2.13P and (ll.2.14p . Thus, in this range and for large A^, effectively the random 
variable Amin in the bipartite problem behaves, in law, as the Wishart minimum 
eigenvalue scaled by a factor A^~^, i.e., Amin — ^ "U^min/A^^- This is also confirmed 
in Eq. (jl.2.5p , where we see that the average trace in the unconstrained Wishart 
ensemble scales as A^^. Thus, in the microcanonical enemble, where the trace 
is constrained to be unity, it amounts to rescale all the Wishart eigenvalues 
by a factor A^^^ for large A^. However, for finite A^, the distributions are very 
different in the constrained and unconstrained ensembles, in particular near the 
upper edge x = 1/N. In other words, the distribution of Amin exhibits strong 
finite size effects. 

1.5 Summary and Conclusion 

In this chapter we have discussed an application of Wishart matrices in an 
entangled random pure state of a bipartite system consisting of two subsystems 
whose Hilbert spaces have dimensions M and A^ respectively with A^ < M. 
The A^ eigenvalues of the reduced density matrix of the smaller subsystem 
are distributed exactly as the eigenvalues of a Wishart matrix, with the only 
difference that the eigenvalues satisfy a global constraint: the trace is fixed to 
be unity. 

We have studied the distribution of the minimum eigenvalue in this fixed- 
trace Wishart ensemble. For the hard edge case (when two subsystems have 
same size M = N), we have shown that the minimum eigenvalue distribution 
can be computed exactly for all A^ in all three interesting physical cases /? = 1, 
/? = 2 and /3 = 4. 



1.5. SUMMARY AND CONCLUSION 21 

What does this exact distribution of Amin tell us about the entanglement en- 
tropy of the bipartite system in a pure state? We have seen before that if Amin is 
close to its maximally allowed value 1/N (set by the unit trace constraint), then 
that configuration is maximally entangled since all the eigenvalues contribute 
equally to the composition of the state. A measure of how close the random 
state is to this maximally entangled state can be estimated by computing the 
net measure of Amin in a small range close to 1/A^, e.g., by the cumulative prob- 
ability Rn{1/N — e) = fiLi_^ Pn{x) dx where e << 1/A^. From our exact cal- 
culation, we see that the pdf of Amin, in all three cases, approaches zero as Amin 
approaches its maximum possible value 1/A^, -P/v(Amin = x,N) ~ (1 — Nxyi^ 
where i^i = {N"^ + N - 4)/2, 1^2 = N"^ - 2 and 1^4 = 2N'^ - N -2. This shows 
that the 'closeness to maximal entropy' measure R]\[{l/N — e) ~ (eiV)'^/'+^. For 
e << 1/A^, this measure is evidently very small. It was argued before |Pag95| , on 
the basis of the computation of the only the first moment of the von Neumann 
entropy (not the full distribution), that a random state is almost maximally 
entangled. Our result shows that the probability that a random state is maxi- 
mally entangled is actually very small. The same conclusion was also deduced 
recently on the basis of the large N computation of the full distribution of the 
Renyi entropy |NadlO) . Thus, the lesson is that conclusions based just on the 
first moment, may sometimes be a bit misleading. 

Here we have discussed only the hard edge M = N case where the two 
subsystems have equal sizes. Our results are of relevance for small systems 
such as when each subsystem consists of identical number of qubits. It would 
be interesting to extend these calculations to the cases when M ^ N. In 
particular, in the context of thermodynamic systems where, for instance, one of 
the subsystems is a heat bath, one needs to study the opposite limit N « M. 
It would be interesting to estimate the distribution of the minimum eigenvalue 
and other measures of entanglement in that limit. 

Finally, we have restricted ourselves here to 'random pure' states where all 
pure states are sampled equally likely. This is the Haar measure. So far, we 
have not discussed dynamics, i.e., the temporal unitary evolution of the sys- 
tem. Under any unitary evolution that is ergodic over the space of all pure 
states, Haar measure is the unique stationary measure of the unitary evolution. 
This ergodicity holds provided one does not have any strict conservation law. 
For instance, the uniform measure over all pure states will not hold under the 
standard microcanonical scenario where the composite system has a fixed total 
energy E (eigenvalue of the Hamiltonian H of the composite system). Only 
those pure states with total energy E, and not all pure states, will be sampled 
by the system under unitary evolution. An appropriate 'stationary' measure 
is then the microcanonical measure which is uniform over all pure states be- 
longing to the fixed E manifold. For such a measure, one can again define the 
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reduced density matrix of the subsystem A and its eigenvalues. It would be 
interesting to study the statistics of the bipartite entanglement (von Neumann 
or the Renyi entropy or the minimum eigenvalue Amin) in this microcanonical 
setting with a fixed total energy E. For such systems, some recent results of 
very general nature (that do not require the detailed knowledge of the Hamil- 
tonian of the system) have been derived |Pop06 fGolOGl iReiOTj which says that 



any pure state, drawn from the uniform measure on the constrained manifold, 
will almost surely be 'maximally entangled' i.e., very close to the maximally 
entangled (centroid) configuration Aj = 1/N for all i, provided M >> N, i.e., 
the environment (subsystem B) is much bigger than the system (subsystem A). 
It would be interesting to compute explicitly the distribution of this 'typicality' 
i.e, the distance between the pure state (drawn from a uniform measure over 
the constrained manifold) and the maximally entangled state for some systems 
with specific Hamiltonians. 
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